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Abstract 

We present the calculation of 0{a 2 s a em ) corrections to heavy-quarkonium total photopro- 
duction cross-sections. Results are given for the colour-octet component of S and P waves. 
The calculation is performed using covariant projectors in dimensional regularization. A 
phenomenological study of the results, including a discussion of the high-energy behaviour of 
the cross sections, is presented. For 7^ energies up to few hundred GeV the NLO corrections 
significantly reduce the scale dependence of the production rates relative to the Born-level 
results. Large small- x corrections arise at higher energies, making the predictions strongly 
dependent on the shape of the gluon density and on the choice of factorization scale. 
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1 Introduction 



The calculation of production and decay processes for heavy quarkonium states has recently been 
put on a solid formal basis by the work of Bodwin, Braaten and Lepage U (BBL). According to 
their results, production and decay rates can be calculated within perturbative QCD as the sum 
of products of short- distance coefficients times long-distance matrix elements. The short- distance 
coefficients are the square of transition matrix elements for production or decay of heavy-quark 
pairs in definite states of colour, spin and angular momentum. The long-distance ones are obtained 
from the matrix elements of quark-antiquark operators with the same quantum numbers as those of 
the short- distance state; these are evaluated between the vacuum and an arbitrary state containing 
the physical quarkonium meson we are interested in, squared and summed over all possible final 
states accompanying the quarkonium. These long-distance matrix elements can in principle be 
calculated on the lattice, and hierarchies among them can be obtained by applying the velocity 
scaling rules of NRQCD pi |3j . Several applications of this formalism have been obtained, and are 
nicely reviewed in ref. [f|, [5| . 

One of the most important consequences of this factorization property of quarkonium produc- 
tion is the prediction that the value of the non-perturbative parameters does not depend on the 
details of the hard process, so that parameters extracted from a given experiment can be used 
in different ones. For simplicity, we will refer to this concept as "universality". Several studies 
of experimental data coming from different kind of reactions have been performed to assess the 
validity of universality. For example calculations of inclusive quarkonia production in e + e~ anni- 
hilation |5[], fixed target experiments 0, 7p collisions || and B decays have been carried 
out within this framework. The overall agreement of theory and data is satisfactory, but there 
are clear indications that large uncertainties are present. The most obvious one is the discrep- 
ancy between HERA data JTT] and the large amount of inelastic J/ip photoproduction predicted 
by applying the colour-octet matrix elements extracted 13, [14| from the Tevatron large-p-r 
data 

In view of this discrepancy, it becomes important to assess to which extent is universality 
applicable. Several potential sources of universality violation are indeed present, both at the 
perturbative and non-perturbative level. On one hand there are potentially large corrections to 
the factorization theorem itself. In the case of charmonium production, for example, the mass 
of the heavy quark is small enough that non-universal power-suppressed corrections can be large. 
Furthermore, some higher-order corrections in the velocity expansion are strongly enhanced at 



the edge of phase-space [T7\. For example, the alternative choices of using as a mass parameter 
for the matrix elements and for the phase-space boundary the mass of a given quarkonium state 
or twice the heavy-quark mass 2m, give rise to a large uncertainty in the production rate near 
threshold. These effects, which are present both in the total cross-section and in the production 
at large-pr via gluon fragmentation, violate universality. This is because the threshold behaviour 
depends on the nature of the hard process under consideration. 

Another source of bias in the use of universality comes from purely perturbative corrections. 
Most of the current predictions for quarkonium production are based on the use of leading- 
order (LO) matrix elements. Possible perturbative i^-factors are therefore absorbed into the 
non-perturbative matrix elements extracted from the comparison of data with theory. Since the 
size of the perturbative corrections varies from one process to the other, an artificial violation of 
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universality is introduced. Examples of the size of these corrections are given by the large impact 
of fc^-kick effects and initial-state multiple-gluon emission in open-charm and charmonium 
production [TP|, |2D . 



In ref. [^TJ we focused on the evaluation of the 0(a 3 ) corrections to quarkonium total hadro- 
production cross-sections. As pointed out in ref. P2| , the impact of NLO corrections can be 
significant and a general study of their effects is necessary. In this paper we concentrate on the 
calculation of the 0(a 2 s a em ) corrections to total photoproduction cross-sections. To carry out our 
calculations, we need a framework for calculating NLO inclusive production cross-sections. As well 
known, the most convenient method for regulating both UV and IR divergences in perturbative 
calculations beyond leading order in a s is dimensional regularization. On the other hand most 
calculations of production cross- sect ions and decay rates for heavy quarkonia have been performed 



using the covariant projection method |23|| , which involves the projection of the QQ pair onto states 
with definite total angular momentum J, and which is specific to four dimensions. 

In ref. we presented a generalization of the method of covariant projection to D = 4 — 2e 
dimensions. In that paper our formalism was shown to provide equal results to calculations 
performed within the "threshold expansion" technique, introduced in |24j, ^5|. In this work we 



apply the covariant-projection technique to the calculation of the 0(aja em ) total cross-sections for 
the photoproduction of several QQ states of phenomenological relevance: , 3 Si and 3 Pj , 
where the right upper index labels the colour configuration of the QQ pair. 

The paper is structured as follows. In Section |2| we briefly review our formalism. A more 
complete discussion can be found in ref. |2T| . Section ^| gives a brief general description of the 



NLO calculation. In particular, we describe the the behaviour of the soft limit of the NLO real 
corrections and the technique used to identify the residues of the IR and collinear singularities and 
to allow their cancellation without the need for a complete D-dimensional calculation of the real- 
emission matrix elements. Section |] presents the various components (real and virtual corrections) 
of the NLO calculation of the production processes. Section |] presents a numerical study of the 
results, with a discussion of the individual components of the cross-sections, of the K factors, 
and of the scale dependence. A discussion of the high-energy behaviour of the production rates 
at NLO and of the uncertainties of the predictions in this regime is included. The last section 
contains our conclusions, as well as a discussion of the relevance of our calculation for the study 
of photoproduction at pr > and z < 1. 

Appendix |A] collects symbols and notations. A summary of all results is provided in Ap- 
pendix [FJ, where the parton-level photoproduction cross-sections are presented in their final form, 
after the cancellation of all singularities. 



2 Introduction to the formalism 

The starting point of the calculation of the production cross-sections is given by the standard 
formula [|TJ: 

da(H + X) = Y J da(QQ[n}+X)(0 H (n)) (1) 
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The quantity (0 H (n)) is proportional to the inclusive transition probability of the perturbative 
quark-pair state QQ[n], with quantum numbers labelled by n, into the quarkonium state H. 
da(QQ[n] + X) is the short- distance cross-section for the production of the perturbative state 
QQ[n]. It can be calculated in perturbative QCD either using threshold-mathing techniques |], [24}] , 
or using projection techniques |E3J. These projection techniques have been recently extended 
for use in .D-dimensions [[H]], to deal with the presence of infrared (IR) and ultraviolet (UV) 
divergences. They will be briefly reviewed here. 

The spin projectors, with non-relativistic normalization for the spinors, for outgoing heavy 
quarks momenta Q = P/2 + q and Q = P/2 — q, are given by 



" f =vib(f-H 7 °(2 + H' (3) 

for spin zero and spin one states respectively. In these relations, P is the momentum of the 
quarkonium state, 2q is the relative momentum between the QQ pair, and m = M/2 is the 
mass of the heavy quark Q. The justification for use of these projectors in D dimensions, and 
a discussion of how to deal with the presence of the 75 matrix, can be found in ref. |2T| . Here 
we limit ourselves to pointing out that the D-dimensional character of space-time is implicit in 
eqs. and appears explicitly when performing the sums over polarizations, as shown later. 

The colour singlet or octet state content of a given state will be projected out by contracting 
the amplitudes with the following operators : 

C\ = —7= for the singlet (4) 



'N, 



c 



C 8 = \f2T^ for the octet (5) 

The projection on a state with orbital angular momentum L is obtained by differentiating L 
times the spin- and colour-projected amplitude with respect to the momentum q of the heavy 
quark in the QQ rest frame, and then setting q to zero. We shall only deal with either L = or 
L = 1 states, for which the amplitudes take the form: 

(6) 
(7) 



As- 


=0,L=0 


= Tr[Ci7 ^l]| 9=0 


Spin singlet S states 




=1,L=0 


= Tr [Cn?A]\ q=0 e a 


Spin triplet S states 


As-. 


=0,L=1 


= ^r{Cn A]\ q=0 ep 


Spin singlet P states 


As-. 


=!,£=! 


= ^-Tr{Cn?A}\ q=0 £ a p 


Spin triplet P states 



(9) 

A being the standard QCD amplitude for production (or decay) of the heavy quark and antiquark 
Q and Q, amputated of the heavy quark spinors. 

The amplitudes As,l wm then have to be squared, summed over the final degrees of freedom 
and averaged over the initial ones. 
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The selection of the appropriate total angular momentum quantum number is done by per- 
forming the proper polarization sum. We define: 

= - 9aP + . (io) 

The sum over polarizations for a 3 Si state, which is still a vector even in D = 4 — 2e dimensions, 
is then given by: 

^e a e* a , = Ii aa , (11) 

Jz 

In the case of 3 Pj states, the three multiplets corresponding to J = 0, 1 and 2 correspond to a 
scalar, an antisymmetric tensor and a symmetric traceless tensor, respectively We shall denote 
their polarization tensors by S^p- The sum over polarizations is then given by: 



£ap££p*> - Jj—^H-apfta'P' (12) 



X^S^S* = 7j[ n o«' n A3' - UaP'^a'p] (13) 
Jz Z 

^£ a p£a'P* = l^-ota'^-PP' + H-aP'Ra'p] ~ Jj _ ^ aP^a'P' (14) 

Jz 

for the 3 Po, 3 -Pi and 3 P2 states respectively. Total contraction of the polarization tensors gives the 
number of polarization degrees of freedom in D dimensions. Therefore 

Nj = E e «4 = n aQ = D - 1 = 3 - 2e (15) 

Jz 

for the 3 S'i state and 

Nj = ( 16 ) 

Jz 

for the 3 Pj states, with 

* = i¥ ^ P-i)(g-2) = (3 _ 2e)(1 _ £)j JVj= (£±^£zS = (5 _ 2e)(1 _ e) . (17) 

The application of this set of rules produces the short-distance cross section coefficients a for 
the ij — > 2S + 1 l/ 1 ' 8 ^ processes: 

da(ij - ^L/> 8] ) = ±Z\^s,l\ 2 d§ , (18) 

s being the partonic centre of mass energy squared. To find the physical cross-sections for the 
observable quarkonium state H these short distance coefficients must be properly related to the 
NRQCD production matrix elements (C^ 8 ]( 2S,+1 Lj)). The cross-sections then read 

- 2S+1 L.r ] - if) = a(i 3 - ^m ^lf^)) > (19) 

^ col^ pol 



4 



where N co i and N po i refer to the number of colours and polarization states of the QQ[ 2 + Lj\ 
pair produced. They are given by 1 for singlet states or N c 2 — 1 for octet states, and by the 
.D- dimensional N/s defined above. Dividing by these colour and polarization degrees of freedom 
in the cross-sections is necessary as we had summed over them in the evaluation of the short 
distance coefficient o. As discussed at length in ref. [jST], our conventions for the normalization 



of the non-perturbative matrix elements differ slightly from the conventional ones introduced by 
BBL0: 

//n \ BBL 

(Oi> = ^ (20) 
(Os) = (0 8 ) BBL . (21) 



3 Soft factorization and the calculation of higher-order 
corrections 

In this section we discuss the role played by the universal IR behaviour of gluon-emission ampli- 
tudes in the calculation of higher-order corrections to total production cross-sections. A consistent 
calculation of higher-order corrections entails the evaluation of the real and virtual emission di- 
agrams, carried out in D dimensions. The UV divergences present in the virtual diagrams are 
removed by the standard renormalization. The IR divergences appearing after the integration 
over the phase space of the emitted parton are cancelled by similar divergences present in the 
virtual corrections. Collinear divergences, finally, are either cancelled by similar divergences in 
the virtual corrections or by factorization into the NLO parton densities. The evaluation of the 
real emission matrix elements in D dimensions is usually particularly complex. In this paper we 
follow an approach already employed in |26|| , whereby the structure of soft and collinear singular- 



ities in D dimensions is extracted by using universal factorization properties of the amplitudes. 
Thanks to these factorization properties, that will be discussed in detail in the following section, 
the residues of IR and collinear poles in D dimensions can be obtained without an explicit calcu- 
lation of the full D-dimensional real matrix elements. They only require, in general, knowledge 
of the D-dimensional Born-level amplitudes, a much simpler task. The isolation of these residues 
allows to carry out the complete cancellations of the relative poles in D dimensions, leaving resid- 
ual finite expressions which can then be evaluated exactly directly in D = 4 dimensions. The 
four-dimensional real matrix elements that we will need can be found in the literature [EL 



3.1 Soft factorization in ggj amplitudes 

At the Born level, the relevant diagrams are shown in fig. [I]. The production amplitude (before 
projection on a specific quarkonium state) can be written as follows: 

A Born = T« D 12 , (22) 

where D\2 = D\ + D2 and the terms D±, D2 correspond to the diagrams appearing in fig. [I] with 
the colour coefficients removed. Using this notation, the amplitude for emission of a soft gluon 
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Figure 1: Diagrams for the 57 Born amplitudes. 



with momentum k and colour label c can be written as follows [p7 |: 



.4, 



soft 



^(T c T a ) 



D 12 +g(T a T c 



'1 j 



ae c Qe c 



ak Qk 



D 12 + 



(23) 



where Q and Q are the momenta of the heavy quarks, and a indicates the momentum and colour 
label of the initial state gluon. Born-level colour-singlet amplitudes vanish, so we will concentrate 
on colour-octet states. Using the projectors defined in the previous section, we can write: 



-4, 



soft 



gy/2D 12 \ Tr (T b T c T a ) 



ae r 



Qe c 

Qk ak 



9^-dJ** 



Qe c Qe c 



+ if 



bca 



Qe c Qe 



ae c Qe c 



ak Qk 



ae c 

Qk Qk ak 



Qk Qk 

where b is the colour label of the colour-octet QQ state. 

In the case of S'-wave production we can set Q = Q = P/2, and we get 

igf bca \Pe 



(24) 



( A m 

^soft 



g=0 



V2 



ae r 



Pk ak 



D 



12 



(25) 



In the case of 3 P-waves we also need the derivative of the decay amplitude with respect to the 
relative momentum of the quark and antiquark. In the soft-gluon limit, we obtain: 



dA. 



[8] 
soft 



dq a 



q=0 



V^2 f . „i, r „ 
9 ~2 I lf 



Pe 



ae r 



Pk ak 



dD 



12 



dq a 



+ d' 



bca 



e a k a 

Pk ~ (Pfcp (PCc) 



D 



12 



(26) 



Choosing a transverse gauge where e c • P = 0, from eq. (9) and the previous two expressions it is 
straightforward to write 



pj 



igf 



bca r 



V2 



Pe c _ aec 
Pk ak 



Tr [CilMBorn] 



g=0 



+ 



gd 



'bca 



V2Pk 



£ a p€^Tr[Cn?ABom]\ q=Q ■ (27) 
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The soft amplitude factorizes if the last term vanishes. This last term can be seen to be propor- 
tional to the amplitude for the production of a 3 Si state, with an effective polarization e e fj given 
in terms of the polarizations of the soft gluon and of the 3 Pj state as follows: 

= &ip(P) ■ (28) 
If the average on the soft gluon D — 1 spatial directions is taken, i.e. 

^££^ = §^MP), (29) 

pol 

one can easily compute the sum over effective polarizations 

/ ^ E = -Nj^n a p(p) , (so) 

where Nj is the number of degrees of freedom of the 3 Pj state in D dimensions. This shows that 
the last term in eq. (0), once squared and averaged on the directions of the outgoing gluon, is 
proportional to the amplitude of a 3 Si state coupled to two gluons , which vanishes by C-parity. 
As a result we obtain the factorized expression: 

2aP M 2 



E l41l 2 = C A g 2 

col, pol 



[ak){Pk) (Pk) 2 



E I^BorJ 2 (31) 
col, pol 



The amplitudes for 'jq — > Qq are IR finite, and there is no need to study their soft behaviour for 
our applications. 



3.2 Kinematics and factorization of soft and collinear singularities 

The kinematics of the process k + p\ — > P + p 2 , where P is the momentum of the heavy quark 
pair, k the momentum of the photon and Pi are the momenta of the massless partons, can be 
described in terms of the standard Mandelstam variables s, t and u: 

s = (fc + pi) 2 , (32) 
t = (P1-P2? = ~(l-x)(l-y) , (33) 

u = (k-p 2 ) 2 = -i(l- x )(l+ y ). (34) 

Here we introduced the Lorentz- invariant dimensionless variables x = M 2 / s and y (— 1 < y < 1), 
defined by the above equations. In the center-of-mass frame of the partonic collisions, the variable 
y becomes the cosine of the scattering angle 9. In terms of x and y the total partonic cross-section 
can be written in D dimensions as follows: 
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a 



JLSJUULf 



'X/X/X/X/N 



ouuuuu 



JUUUU 




b 



c 



Figure 2: Diagrams for the real corrections to the 57 channels. Permutations of outgoing 
gluons and /or reversal of fermion lines are always implied. 



where Ai = is the spin- and colour-averaged matrix element squared in D dimensions and 

d$( 2 )(x, y) is the D-dimensional two-body phase space: 

A e /A 7r \ e 1 

d$ {2) (x,y) = - [-) T(l + e) — (l-x)^(l-y 2 rdy, (36) 

with 



,7T 2 



K = T(1 + e) r(l - e) = 1 + e 2 — + 0(e 3 ) . (37) 

o 

The soft and collinear singularities are associated to the vanishing of t or u, which appear at most 
as single poles in the expression of M.. One can therefore introduce the finite, rescaled amplitude 
squared M: 

M = ^= sm-xni-y*) ™- < 38 » 

In terms of A4, the partonic cross-sections read as follows: 

4(7 r 1 
a(x) = / dy(l-y 2 )- 1 -*M(x,y)dy, (39) 

S J — 1 

Soft and collinear singularities are now all contained in the universal poles which develop as x — > 1 
and y 2 — > 1. The residues of these poles can be derived without an explicit calculation of the 
matrix elements, as they only depend on the universal structure of collinear and soft singularities. 
We will carry out an explicit evaluation of these residues in the next section. 



4 Results 

4.1 91 -> gO® processes 

We start by considering the soft limit, x — > 1. The following distributional identity holds for small 

a-,,- - J^ (1 -, ) + (^) p - 2£ (^) /0( ,) (41) 
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where p = M 2 / S-yh, h being the initial state target hadron, and (3 = a/1 — P- The p-distributions 
are defined by: 

J dx [d(x)) p t(x) = J dxd(x) [t(x) - t(l)] . (42) 
We can therefore write, with obvious notation, 

a( x ) = a x=1 + a x+1 . (43) 
The first term on the right-hand side is given by the following expression: 

AC r~ 4e r 1 

a x=1 = -— — 6{l-x) j_ i dy(l-y 2 )- 1 ^M(x = l,y). (44) 

The x — > 1 limit of M. can be easily derived from eq. (PT|): 

M(x,y) ^ sg 2 C A (l-y) 2 M Born , (45) 

where A^Bom is the D- dimensional Born amplitude squared for the — > process, which 
is independent of y. The integration over y of eq. (|44|) is elementary, and leads to the following 
result: 

(T x= l = \ I ^2 HC A — CTBorn , (46) 

where a s is related to the D-dimensional bare coupling af 1 and to the renormalization scale p by 
a s p 2e = = g 2 /An. H is defined by 

g = r(i + iy-20 = 1 + t + ^-T^ + o <^ (47) 

and UBorn is the .D-dimensional Born cross-section: 

(TBorn = 71 / M 4 M B om5(l ~ x) = <T 6(1 - x) . (48) 

The collinear singularities remaining in <j x -l\ can be factored out by using the following distribu- 
tional identity: 



where the distributions on the right-hand side are defined by: 



+ 0(t) , (49) 



The contribution a x ^\ can then be split into two terms: 

(J X -L1 = <J y =X + CJfinite • (51) 
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The term cr n mte has no residual divergences, and is given by the following expression: 



2C 



O finite 



1 



1 — X 



-1 



dy 



1 



1 



1-3/, 



1 + 1/, 



M{x,y) . 



(52) 



We removed the + from the 1/(1+2/)+ distribution because no collinear singularity can arise from 
the emission of the gluon collinear to the photon, and A4(x, —1) = 0. <7fi n i tc explicitly depends on 
the nature of the quarkonium state produced. For processes whose Born contribution vanishes, 
this is the only non zero term. u y =i is given by: 



4C 4" 



2e 



1 — x 



2e 



loefl 



1 — X 



M{x,y= 1) . 



(53) 



The limit for y — > 1 of A4(x,y) is universal, thanks to the factorization of collinear singularities: 

y-,1 _ ^ „ , s 1 — X 



M(x,y) 



ITT Sa^Pgg[X) 



-M 



x 



Born 



Using this relation we get: 



si ^ P 99( X ) i 1 ~ X ) X a 



1 — X 



2e 



log(l — x) 

1 — X 



with 



and 



x 



1 — X 



2C A —— + z — - + x(l-x) 
.1 — X X 

1 1 

r = - - 7e + log(47r) 



(54) 



(55) 



(56) 



(57) 



The collinear poles take the form dictated by the factorization theorem. According to this the 
partonic cross-section can be written as: 



da 7j (p^,p h ) = ^d(T 7/ (p 7 ,xp /l )r ij (x)d 



x 



Tij(x) = 8ij8(l — x) 



e 2tt I ^| 



P 4i (x) + /x, ,(./•) 



(58) 
(59) 



where da is free of collinear singularities as e — ► 0. Here we allowed the factorization scale /xp 
to differ from the renormalization scale fi. The functions Vij(x) are the D = 4 Altarelli-Parisi 
splitting kernels, collected in Appendix [A], and the factors Kij are arbitrary functions, defining 
the factorization scheme. In this paper we adopt the MS factorization, in which Kij(x) = for 
all For the definition of Kij(x) in the DIS scheme, see for example ref. [2S]. 

(c) 



Expanding eq. (|58f) order-by-order in a s , we extract the counter-term (Ty=i, defined by: 



a 



(c) 
y=l 



2 \ € 



(60) 
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'X/X/V/V/^r 



a 




b 



Figure 3: Diagrams for the *yq channels. Reversal of fermion lines is always implied. 



with 



2C A 



x 



(1-x), 



+ 



X 



+ x(l — x) 



X 



+ {b + 4C A log/3) 5(1 -x). 



(61) 



Putting all pieces together, we come to the final result for the real-emission cross-section: 



Z7T 



+ 



n 2 \ 2 

y-41ogr + 81og 2 rj - ^n f T F 



xV gg (x) log 4 + 2^(1 - x)P gg (x) f 10 ^ 1 X > 
/ip y 1 x 



+ 



x 



6{l-x)+ 

J\ 9 [Q [8] ](x) 



\Q [> 



1q[8] 3p[8] 3p[S 
J -M) > *2 



(62) 



where f e (s) is defined in Appendix [A] and crjf [7<? — > <2' 8 '] is the D-dimensional, Born-level partonic 
cross-section for the production of the quarkonium state H via the intermediate state, after 
removal of the 5(1 — x) term (see eq. (^8|)). The finite functions f ig (x), obtained from the explicit 
evaluation of eq. (|52|), are collected in Appendix [B|. 



processes 

The Born-level processes <27 — ► identically vanish. As a result IR divergences at 0(a 2 a em ) 
and virtual corrections are not present. The only possible singularities appearing at this order 
come from the emission of the final-state quark collinear to the initial-state one or to the photon.. 
The behaviour of the amplitude (fig. 0(a)) for Q' 8 ' = 1 Sq'\ 3 Pq 8 \ 3 P^ m the y — > 1 collinear limit 
is again controlled by the Altarelli-Parisi splitting functions: 

M(x, y) y -^87rs a ^ P gq (x)^—^M Bo m • (63) 
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In analogy to the 75 case, one introduces the following counter-term in the MS scheme: 



a 



(c) 



1 = (tt) P OT (a;)a:«7o 

e 2vr V^f/ 



(64) 



where V gq (x )is defined in Appendix [A| 

Following a procedure analogous to the one detailed in the case of production, we find the 
following result: 

a, 



a H [ iq ^Q [8] q](x) 



71 



Ml 



-P 9<? (x)loj 



s(l — 2) 



.r 



[fi 



[8] _ 1 q[8] 3 p[8] 3 p [8 



(65) 



where the functions f iq {x) are collected in Appendix ||, together with the result for 3 Pi production, 
for which no collinear singularity is present to start with. 

For 3 S'{ 8 ' production only the diagram in fig. ^(b) gives a non vanishing contribution. In the 
y — > 1 collinear limit we have: 



M(x,y) 2=3 &Ksa®($ J P rr {x)- — -M Born . 



1-x 



X 



(66) 



where A^Bom is the Born amplitude for the process qq — > 3 Sj 8 ' and -P ?7 (x) is the Altarelli-Parisi 
photon-splitting function in D dimensions. Introducing the following counter-term in the MS 
scheme: 



1 CVem&f 



A7) V q ry(x)xa 



2/=l o 12 

where V qi [x) is defined in Appendix 0, we can write the following partonic cross-section: 
a H [ iq ^ 3 S?q](x) = ^La»iqq^ 3 S [ * h 



(67) 



7T 



X 



-P g7 (x)log " - 



L> F x 2 (l -x) 



+ / 7? m(^) (68) 



Having absorbed the collinear 7 — > g divergence in the photon structure function, the total cross- 
section will now include a piece proportional to the resolved component of the photon (see, e.g. , 
ref. J2g|T§): 



da jh (pi,p h ) = ^2da y i(xS,iiF,^)Fi h (x,n F )d 



x 



+ d<T fc/ (xip 7 , x 2 Ph, /i F )F fc7 (xi, fJL 1 )Fi h {x 2 , n F )dx 1 dx 2 
k,i 



(69) 



where -Ffc 7 (x, /i 7 ) is the density of the parton k in the photon, and Fkh(x, fip) is the density of the 
parton k in the hadron h. The full 0(ct 3 ) expression for the parton-parton 3 5^ cross-sections to 



be used in the resolved-photon contribution to the above equation can be found in [21 
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Figure 4: Feynman diagrams contributing to the one-loop corrections to the processes 
57 -> QQ[ 3 P,S 81 ], 57 -» QQP-fj 81 ] and 57 - QQpS? 1 ]. 
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Table 1: "Diagrammatic" partial virtual QCD corrections to the processes — > [ 1 Sq ] 
(diagram multiplicities are included) 
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Table 2: "Diagrammatic" partial virtual QCD corrections to the processes 57 
(diagram multiplicities are included). 
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Table 3: "Diagrammatic" partial virtual QCD corrections to the processes — ► [ 3 P^ ] 
(diagram multiplicities are included). 

4.3 Virtual corrections 

We present in this section the results of the calculation of the 1-loop diagrams necessary for 
the evaluation of the virtual corrections to the production matrix elements. These results can 
be obtained in a straightforward way from the calculation of the virtual corrections to 1-loop 
hadroproduction, presented in ref. We shall therefore limit ourselves to presenting the final 
answers. The relevant Feynman diagrams are shown in figures |], and the results are given 
diagram by diagram in tables HI 0, [| In these tables we report the contribution of each diagram 
k, indicating separately the colour factors f k . The expressions Vk appearing in the tables are 
defined by the following equation: 

a H {ij ^ Q] = [y _ Q)f e ( s )J2V k f k 5(l-x). (70) 

71 k 

where the sum extends over the set of diagrams and / e (s) is defined in Appendix |A]. 

The singularity structure of the virtual corrections is dictated by the renormalization properties 
of the theory, by the universal form of the Coulomb limit, and by requirement that soft and 
collinear singularities cancel against the real corrections evaluated above. The form of the virtual 
corrections to the cross-section is therefore the following: 

^ = a ^f e (s) x (A. + {Cf - ±C A f- -cj + + -^n f T F + 2D®\, (71) 

2tt [euv v \ef R Qe m J 3e m J 

where we explicitly labelled the e's to indicate their origin, and where all of the state dependence 



15 



is included in the finite factor Dq. The quark-antiquark relative velocity 2v and 60 are defined in 
Appendix |A|. rif is the number of flavours lighter than the heavy, bound one. 

Summing the contribution of all diagrams, we obtain the following results for the colour-octet 
coefficients Dq : 

< = *{-*+t) +c 4 + $ (72) 

< = -4C F + C. 4 (f + ^ + £). (74) 
The final results for the finite sums of real plus virtual corrections are collected in Appendix [B]. 



5 Phenomenology 

In this section we study some of the properties of the higher-order corrections calculated in this 
paper, and their effect on typical Born-level predictions. A more thorough phenomenological study 
including a comparison to currently available data and fits to the non-perturbative parameters 
will be presented elsewhere. 

To start with, we show in fig. |5| the comparison between Born and NLO total photoproduc- 
tion cross-sections as a function of the photon beam energy, in the energy range of fixed-target 
experiments. Here and in the following, we concentrate on the colour-octet contributions of 1 Sq, 
3 Pq and 3 P2 states, which are by far the dominant processes. We fix m c = 1.5 GeV and we take, 
for the sake of definiteness §, (0\Oi{ x S )\0) = 0.01 GeV" 3 and (0\Ot( 3 P )\0)/m 2 c = 0.01 GeV" 5 , 
with (0\OS(*Pj)\0) = (2J+ 1)(0|^'( 3 P )|0). 

The scale dependence of the results is displayed by the curves relative to the scale choices 
A*r = A*f = °>M, with a = 1/2, 1 and 2. We chose the MRSA set of parton densities |JD|], including 



the low-Q 2 corrections |[T| necessary to evaluate consistently the cross-sections when using the 
smaller choice of scales, fip = M/2. These low-Q 2 corrections are essential to properly estimate 
the true scale-dependence of the calculation. Similar curves for the energy range typical of the 
HERA collider (given as a function of S 1P , the 7p CM energy), are shown in fig. Notice the 
significant improvement in the scale dependence when the NLO corrections are included. This 
improvement is less remarkable at low energies. The reason for this is that the 57 — ► qQ is 
particularly important at low energies. This process first appears at 0(a 2 a em ), and is therefore 
calculated at LO only. Its contribution becomes less important at high energies, where the gluon- 
initiated process dominates, and the 0(a 2 a em ) calculations are therefore genuinely NLO. As a 
result of this, the if-factor (defined as the ratio of the 0(a 2 a em )/Born results) is larger at low 
energy (exceeding a factor of 2) than at high energies. 

The relative importance of the gluon and quark processes is displayed in figs. [7] and || for 
the fixed-target and collider-energy ranges, respectively. Notice that the contribution of the qj 
channel is always negative. This is the result of the subtraction of the mass singularities present 
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Figure 5: Born and NLO total photoproduction cross-sections as a function of the 
photon beam energy. 
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Fi gure 7: Separate contributions from Born, 0(a s a em ) *yg and 0(a s a em ) *yq processes 
to total photoproduction cross-sections as a function of the photon beam energy. 



in the q^y — > qQ channel when the final-state quark is emitted collinear to the beam. The collinear 
singularities are absorbed in the gluon parton density at NLO, as explained in the previous sections. 
The net result of this subtraction, in the MS scheme, is a negative contribution of the factorized 
(77 process. 

An important element in previous Born-level extractions @, |9|] of the non-perturbative colour- 
octet parameters from the fit to fixed-target and HERA data is the set of relations: 

^p pM] = 3 mg<0|Og(3^)|0> M M = mg<0|Ogf(3n)|0) M 

1 J (OlOtCSo)^} 1 2 J (OlO^^o)^) 1 J v ; 

In these relations we assumed, as usual, (0|C^( 3 Pj)|0) = (2 J + 1)(0|C^( 3 P )|0). As a result of 
these relations, one could use the following identity: 

aW=0a[ 1 4 81 ], with e = (0|^(%)|0) + 7 (0| ^ (3 2 jPo)l ° ) . (76) 

The validity of these Born-level relations when the NLO corrections are included is studied in figs || 
and [II], which show the ratios of 3 Pj 8 ^ and 1 5q 8 ' production at NLO. As the figures indicate, the 
relation in eq. (|76| ) holds at NLO to within 10%. Notice that the individual relative contributions 
of the 3 P(S 81 and 3 Pj 8] states can however change by up to 30% with respect to the Born-level 
prediction a[ 3 P^ ] ]/a[ 3 P?] = 4/3. 
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Figure 8: Separate contributions from Born, 0(a 2 a cm ) 75 and 0(a 2 a cm ) 7g processes 
to total photoproduction cross-sections as a function of the 7p CM energy. 



5.1 High-energy behaviour 

While the current experiments only allow to study ip photoproduction up to CM energies of the 
order of few hundred GeV, it is interesting to consider the behaviour of the total cross-sections in 
the asymptotic regime. In this regime, interesting phenomena are expected to take place, because 
of the potentially large small-x effects associated to the presence of diagrams with t-channel gluon 
exchange. Because of these contributions, the total partonic cross-sections tend to a constant limit 
when x — > 0. It is easy to see, in fact, that: 



a% i -> Qi\ -(-) Z, Z Q a»[ 19 -> Q] , (i = q,g), (77) 



7T 



where 
and 



Z q = Cf , Z g = Ca , (78) 
A/r2 

(79) 
(80) 
(81) 
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Figure 9: Ratios of the 3 Pj ' to 1 S\ j ' production rates at NLO, normalized to a common 
choice of non-perturbative matrix elements, as a function of the photon beam energy. The 
Born-level predictions for these ratios are equal to 3 and 4 for the ratios of 3 Pq 8 ^ / 1 Sq'^ 
and 3 P^/ 1 S'q 8 ', respectively. 
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Figure 10: Same as fig. 0, as a function of the 7p CM energy. 



Notice that in the x — > limit the ratio between the quark and gluon channels is equal to Cf/Ca, 
as expected. The factors for the 'jq channel are the same as those for the gq channel ||21|| . The 



factors for the gammag channel are half those for the gg channel [pi]] . This is as expected, since 
in this last case the soft t-channel gluon can be radiated from either of the two initial state gluons. 
These results are trivial but useful cross-checks of our calculations. 

Notice also that the small- a; partonic cross-sections tend to a negative value, unless the factor- 
ization scale is chosen to be very small. The large contribution coming from the small-x region, 
and the large scale-dependence of the asymptotic x — > limit, suggest that the NLO calculations 
should display a large dependence on the scale and on the shape of gluon densities at sufficiently 
high CM energies. A similar behaviour has already been observed in the NLO hadropro duct ion 
case [ p2]| . 

To study this issue in more detail, we assume the gluon density to take, at a given scale /xp, 
the form: 

G(x) = , < 5 < 1 . (82) 



x 



1+5 



The softer the gluon density (5 — ► 0), the more important the small-x contributions will be. In 
the extreme case of 5 = 0, it is easy to find the following result for the NLO/LO ii'-factors of the 
total production cross setions, exact up to order p = M 2 / S lp : 



K[Q] 



1 + — 

7T 



1 11 M 2 

-A tot [Q}+C A k Q --C A \og{—\ 

Z jjp 



- C A Z Q log(-) 



(83) 



where the coefficients A to t[Q] are collected in Appendix IB], the -Zq's are given in eqs. (|79])-(pl 
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Figure 11: High-energy behaviour of the K factor for the 75 process, using a gluon 
density G{x) = 1/x, function of the 7p CM energy. 
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As anticipated, large negative logarithmic terms arise. The scale dependence of the coefficient 
of the log p terms cannot be used to change the overall sign, since it is formally compensated by 
the scale dependence of the gluon density, which in the 5 = case is given by: 



dG(p) 
d log 



2^ 



1 dz 



Pgg(z) G(p/Z,H F ) 



Ca^s 1 

G(p) log-. 

7T p 



17) 



The behaviour of these functions is shown in fig. |TT| The fT-factor becomes negative already for 
relatively small values of JS JP , making the NLO estimates unreliable. 

This result is not inconsistent with the nice behaviour of the cross-sections at \ / S 1P ~ 



100 GeV 

seen in fig. |], since in that figure the results were obtained using the more realistic set of parton 
densities MRSA. Typical values of 5 at /if = 3 GeV for the current fits of the gluon densities |30]. [52 
are around 5 



0.3 (see fig. [12]). The stability of the predictions in fact improves for 5 > 0, 
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Figure 12: Small- x behaviour of the gluon density at /ip = 3 GeV, for the MRSA 



and CTEQ4M [32| fits. Arbitrary normalization. 



when the contribution of the small-x regime is suppressed. One can easily verify in fact that for 
<51og(l/p) ^> 1 the asymptotic behaviour of the i\"-factor is given by: 

K[Q] p ^ l + ^- (Z Q + O(S,p)) , (88) 

7T 

Depending on the value of 5, the NLO correction can still be large and negative. Figure ^ shows 
the high-energy behaviour of the A"-factor for the values 5 = 0.2 and 5 = 0.4 Q. As predicted by 
eq. (^), the A"-factor tends to a constant. 



6 Conclusions 

We presented in this paper the first calculation of 0(a 2 s a em ) corrections to the total quarkonium 
photoproduction cross-sections. The calculations have been performed in the framework of the 
NRQCD approach to quarkonium production. As in the Born-level case, the production is domi- 
nated by colour-octet QQ states. The contribution to the cross-sections from the NLO corrections 
is large at typical fixed-target energies (up to a factor of 2 increase over the Born results), and 
decreases at energies typical of the HERA collider. In this energy range the NLO corrections 
significantly improve the stability of the calculated rates under variations of the renormalization 

4 We chose not to use directly current fits of the gluon densities, since in any case these are not reliable in the 
range x < 1CT 5 (y/S^k, 1 TeV). 
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Figure 13: High-energy behaviour of the K factor for the 75 process, using a gluon 
density G(x) = l/x 1+s , as a function of the 7p CM energy and for different values of 5. 



and factorization scales. For energies above few hundred GeV in the 7^ CM frame, large and 
negative small-x contributions dominate the production rate, and make the NLO evaluation of 
the total cross-sections strongly dependend on the scale and shape of the gluon density, calling 
for the resummation of small-x effects. 

The calculations presented in this paper do not directly improve our knowledge of the quarko- 
nium photoproduction at p T > 0, since production at pt > is a process of 0{a 2 s a em ) at the Born 
level. Nevertheless they provide a fundamental element in the determination of the range of px 
and z = p^-Pn/p^-Pn where the Born calculations can be reliably used. In fact the large negative 
contributions arising in the pr = and z = 1 regions from the virtual corrections evaluated in 
this paper affect, via perturbative Sudakov effects, the estimate of the production rates near the 
end-points of the elastic region. The implications of these effects in the comparison of the z and px 
distributions measured at HERA with the predictions of QCD will be examined in a forthcoming 
publication. 
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A Symbols and notations 



This Appendix collects the meaning of various symbols which are used throughout the paper. 
Kinematical factors: 



M 2 M 2 , 

M = 2m, v = \ l , p= — , p = {l-p)* , 



(89) 



"iv 



where s is the partonic center of mass energy squared and S 7P is the 7-hadron one. v is the velocity 
of the bound (anti)quark in the quarkonium rest frame, 2v being then the relative velocity of the 
quark and the antiquark. We also define: 



fe(Q 2 ) = (^Pj r(l + e) = 1 + 6 i- lE + log(4vr) + log £\ + 0(e 2 ) , 



(90) 



and we denote a perturbative QQ state with generic spin and angular momentum quantum num- 
bers and in a colour-singlet or colour-octet state by the symbol 



= QQ[^L/' S] ] 



(91) 



Altarelli-Parisi splitting functions. Several functions related to the AP splitting kernels enter in 
our calculations. We collect here our definitions: 



where 



D F 



V qi (x) 
P 9g (x) = C F 

V gq (x) = C F 



x 2 + (l 

2 
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x z + (1 - x f 

l + (l-a:) 2 
x 

1 + (1 -xf 



.1 r - 2e x(l — x) 



e x 



x 



PggW 



1C A 



x 



+ 



x 



+ x(l — x) 



P 99 (x) = 2C A 



■ X 
X 



X 



(1-x), 



1 — X /-, X 

H h xll — x) 

X 



V 99 (x) = p 99 {x) + ( b o + 4C A log (3) 5(1 - x) 



11 2 
bo = yC A - -T F n f 



(92) 
(93) 

(94) 

(95) 

(96) 

(97) 
(98) 

(99) 



with rif the number of flavours lighter than the bound one. The Pij are the D-dimensional splitting 
functions which appear in the factorization of collinear singularities from real emission, while the 
functions Vij are the four- dimensional AP kernels, which enter in the MS collinear counter-terms. 
The p-distributions are defined by: 



y 1 dx [d(x)] p t(x) = ^ dxd(x) [t(x) - t(l)] 



(100) 
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Colour Algebra 

[T a ,T b ] = if abc T c (101) 

tab 

{T\T b } = d abc T c + — (102) 

Ti(T a T b ) = T F 5 ab T F = \ 

J2a(T a T a ) i:j = C F S i:j C F = = | 

E bc f abc f ebc = C A S ae C A = N C = 3 

E abc d abc (T a T b T% = C 2 {F)5 l0 C 2 {F) = Egzgg£z3 = f 

D F = £<f« = JV c = 3 (104) 

i 

D A = ]T 5 aa = N 2 C - 1 = 8 (105) 



The following formulas were found to be useful: 



£l£Z£ = Ifak-j^Su) (106) 

Tr(T a T b T c ) = h,d abc + if abc ) (107) 

Tr(T a {T b ,T c }) = ^d abc (108) 

CpDp = T F D A (109) 



(Oi) 
(Os) 



Notice that, according to the discussion in ref. [21 
introduced in ref. [p]] (and labelled here as BBL): 



our conventions differ slightly from those 



(Oi) 



BBL 



2N r 



(Os) 



BBL 



(110) 

(111) 



A.l 0(a s a em ) cross-sections 

The D-dimensional cross-sections read 

- ^LjM - if) = a(ij - ^L™) ^"f^^ , (112) 

' coi-' ' pol 

the short distance coefficients a having been calculated according to the rules of Section ||. N co i 
and Np i refer to the number of colours and polarization states of the QQ[ 2S+1 L^' 8 ^] pair produced. 
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They are given by 1 for singlet states or Da = N 2 — 1 for octet states, and by the D-dimensional 
Nj's defined in Section ||. Recall that the matrix elements appearing in the equations above are 
meant to be the bare D-dimensional ones. Making use of their correct mass- dimension, 3 — 2e 
and 5 — 2e for S and P wave states respectively, gives the right dimensionality to .D-dimensional 
cross-sections, i.e. 2 — D = — 2 + 2e. 

We shall use the short-hand notation 

a H (ij -> Q^) = a(tj -> -> H) (113) 



to indicate the production process of the physical quarkonium state H via the intermediate QQ 
state Q^ = QQ[ 2S+1 Lj m }. 

The D-dimensional Born cross-sections read: 



<ra(97 _ = ^4^ 1-2. _ (^o) (114) 

*.<rr - ^^(i-.U) ^ - ^ duo 

<™(97 -» 3 -Pi [S| ) = (117) 
BornVy ' ; m 7 (l-e)(3-2e) V ; D A (1 - e)(5 - 2e) V ; 

B Summary of 0(a 2 s a em ) Results 

We define: 

^{13 - Q)6(l -x) = < rn ( 7 j - Q) (119) 
The 0(a 2 a em ) cross-sections are given as a function of the variable x = M 2 /s. 
The 7# — > Q} 8 *X channels 



a « bg ^s S ? g ] = ±^- !lg fsf](x)(Of(^)) (120) 



+ (10x 2 + 4x 3 + 2x 4 ) logx 
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3(-l + :e) 4 (1 + :e) 



5 - 4x - 71x 2 + 26x 3 - 67a; 4 + 75x 5 



121^ 



^bg^P^g] = ^^ / 7g N 8] ](x)(^( 3 Pr)) (122) 
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+247a; 6 - 23a; 7 - 82a; 8 - 71x 9 - 32a; 10 - 3x n + (-18a; 2 - 6a; 3 
-270a; 4 + 138a; 5 - 6a; 6 + 186x 7 + 186a; 8 + 66a; 9 + 12a; 10 ) logx] (123) 



Q [8] 9] = <N - Q [8] ] (5(1 ~ x) + g {A tot [QM] 5(1 
xP flfl (x) log 3- + 2x(l - a;)P ra (a; 



x) 



+ 



+ 



2tt 

4m 2 „ s (\og{l — x) 



1 — x 



X/jL f 



1 — X 



+ 



(124) 



c,(-io + £) + c(»-£) 

-AC A log /3 + 8C A log 2 [3 + 2 b log + 8C A log /3 log — 

/ip A*F 



C F - 



14 



7T 

T 



C„ 3- 



7T 

IF 



2m 



-AC A log /3 + 8C A log" + 2 b log ^- + 8C A log /3 log — 
-8CV +C A ^ + ^ + lo g 2 N ) 



-4C A log + 8C A log 2 /3 + 2 6 log 



lC A logp\og 



2m 



2C A 



-l-x 1 -2x 6 + 2x b + x b + x* 



(1 -x)(l + x) 3 

+ (1 - 4a; 2 - 2a; 4 - Ax 6 + x s ) logx 

— - [-43 + 14x + 193a; 2 - 82a; 3 - 457a; 4 + 342a; 5 

27(1 - xf(l + xf L 

+124x 6 - 302a; 7 + 113a; 8 + 28a; 9 + 43a; 10 + 27a; 12 + (27 - 90a; 2 + 42a; 3 
+135x 4 + 78a; 5 - 750a; 6 + 270a; 7 - 117a; 8 - 6a; 9 - 192a; 10 + 27a; 12 ) log a; 



18(l-a;) 3 (l + x) E 



-53 + 16a; + 119a; 2 - 86a; 3 - 521a; 4 + 621a; 5 - 55a; 6 



(125) 



(126) 



(127) 



(128) 



-97a;' + 394a; 8 - 445a; 9 + 80a; 



10 



9a; 



ii 



36a; 12 + (36 - 126a;" + 663;^ 



+210x 5 - 798a; 6 + 1242a; 7 - 810a; 8 + 402a; 9 - 168a; 10 + 36a; 12 ) log a; 



(129) 

378a; 4 
(130) 
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The 'yq(q) — > Q' 8 'X channels 



7T 



X 



x r> ( M 4m2 (!- 



x 



+ L> F x 2 (l - x) 



+ AJ 3 M 8] ](*)|> (131) 



2 2 2 



(2m)' 



(132) 



where 



A 9 N 8] ](*) 

Aat^o ^]( x ) 



71 



X 



Pgq(x) lof 



4m (1 — x) 

2 

X jjj-p 



+ C f 



+ f, q [Q [8] )(x) , 



[e M = i 5o M S P M 3 P M] 



4 

1 128 



x(-l + x)(l + 3x) 
1 



(— 1 + x)(— 5 + 4x + 4x ) — x logx 



C F (-l + x)(l - logx) 



9 

Cf 
12 



-(-1 + x) (43 - 14x + 4x 2 ) + (9 - 9x + 4x 2 ) logx 
o 

-(-1 + x)(53 - 16x + 20x 2 ) + (12 - 12x + 5x 2 ) logx 
3 



(133) 



(134) 
(135) 
(136) 
(137) 

(138) 
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